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We analyze the effect of plasma screening on nuclear reaction rates in dense matter composed of 
atomic nuclei of one or two types. We perform semiclassical calculations of the Coulomb barrier 
penetrability taking into account a radial mean field potential of plasma ions. The mean field po- 
tential is extracted from the results of extensive Monte Carlo calculations of radial pair distribution 
functions of ions in binary ionic mixtures. We calculate the reaction rates in a wide range of plasma 
parameters and approximate these rates by an analytical expression that is expected to be applica- 
ble for multicomponent ions mixtures. Also, we analyze Gamow-peak energies of reacting ions in 
various nuclear burning regimes. For illustration, we study nuclear burning in 12 C- lfi O mixtures. 



I. INTRODUCTION 

Nuclear fusion in dense stellar matter is most impor- 
tant for the evolution of all stars. The main source of 
energy for main-sequence stars is provided by hydrogen 
and helium burning. Subsequent burning of carbon and 
heavier elements [1[ drives an ordinary star along the 
giant/red-giant branch to its final moments as a nor- 
mal star. Nuclear burning is also important for compact 
stars. It is the explosive burning of carbon and other 
elements in the cores of massive white dwarfs; it triggers 
type la supernova explosions (see, e.g., Q and references 
therein). Explosive burning in the envelopes of neutron 
stars can produce type I X-ray bursts [H and superbursts 
(that are observed from some X-ray bursters; e.g., Refs. 
0) HI)- Accreted matter which penetrates into deeper 
layers of the neutron star crust undergos pycnonuclear 
reactions. They can power thermal radiation observed 
from neutron stars in soft X-ray transients in quiescent 
states (see, e.g., Refs. BiSilllSMIllIll)- Thus, 
nuclear reactions are important in stars at all evolution- 
ary stages. 

It is well known that nuclear reaction rates in dense 
matter are determined by astrophysical S'-factors, which 
characterize nuclear interaction of fusing atomic nuclei, 
and by Coulomb barrier penetration preceding the nu- 
clear interaction. We will mostly focus on the Coulomb 
barrier penetration problem. Fusion reactions in ordi- 
nary stars proceed in the so called classical thermonu- 
clear regime in which ions (atomic nuclei) constitute a 
nearly ideal Boltzmann gas. In this case the Coulomb 
barrier between reacting nuclei is almost unaffected by 
plasma screening effects produced by neighboring plasma 
particles. The Coulomb barrier penetrability is then well 
defined. 

However, in dense matter of white dwarf cores and 
neutron star envelopes the ions form a strongly non-ideal 
Coulomb plasma, where the plasma screening effects can 
be very strong. The plasma screening greatly influences 
the barrier penetrability and the reaction rates. Depend- 
ing on density and temperature of the matter, nuclear 



burning can proceed in four other regimes [14| . They are 
the thermonuclear regime with strong plasma screening, 
the intermediate thermo-pycnonuclear regime, the ther- 
mally enhanced pycnonuclear regime, and the pycnonu- 
clear zero-temperature regime. The reaction regimes will 
be briefly discussed in Sec. Ill Al In these four regimes the 
calculation of the Coulomb barrier penetration is com- 
plicated. There have been many attempts to solve this 
problem using several techniques but the exact solution 
is still a subject of debates (see [HI, [HI, QjJ and reference 
therein). 

In our previous paper (l7j we studied nuclear reactions 
in a one component plasma (OCP) of ions. Now we ex- 
tend our consideration to binary ionic mixtures (BIMs). 
In Sec. |TT] we discuss physical conditions and reaction 
regimes. In Sec. Mil we describe the results of our most 
extensive and accurate Monte Carlo calculations, ana- 
lyze them and use to parameterize the mean-field poten- 
tial. In Sec. IIVI this potential is employed to calculate 
the enhancement factors of nuclear reaction rates. We 
study the Gamow peak energies and the effects of plasma 
screening on astrophysical S'-factors in Sec. |V] Section 
I VII is devoted to the analysis of the results. We conclude 
in Sec. I VIII In the Appendix we suggest a simple gener- 
alization of our results to the cases of weak and moderate 
Coulomb coupling of ions. 



II. PHYSICAL CONDITIONS AND REACTION 
REGIMES 

A. Nuclear reaction regimes 

We consider nuclear reactions in dense matter of white 
dwarf cores and outer envelopes of neutron stars. This 
matter contains atomic nuclei (ions, fully ionized by elec- 
tron pressure) and strongly degenerate electrons, which 
form an almost uniform background of negative charge. 
For simplicity, we will not consider the inner crust of 
neutron stars (with density higher than the neutron drip 
density ~4x 10 11 gem -3 [lg|), which contains also free 



2 



degenerate neutrons. 9.5 

We will generally study a multicomponcnt mixture of 
ions j = 1,2,... with atomic mass numbers Aj and 
charge numbers Zj. The total ion number density can 9.0 
be calculated as n; = y\ rij, where rij is a number den- 
sity of ions j. It is useful to introduce the fractional r- — ■ 



number Xj — rijjn\ of ions j. Let us also define the av- fcEj . 8.5 - 



erage charge number (.Z) 



y\. Xj^j and mass number 



(A) = yj. XjAj of ions. The charge neutrality implies 
that the electron number density is n e = (Z)n{. The elec- 
tron contribution to the mass density p can be neglected, 
so that p w (A)m u n, where m u w 1.66054 x 10~ 24 g is 
the atomic mass unit. 

The importance of the Coulomb interaction for ions j 
can be described by the coupling parameter Tj 



o 

o 




7.5 



r, 



Z/ a P . 



1/3 



7.0 



Carbon ignition — ^ x 

Tij = 1 year x 



r 



>j 



10 6 year , A X 
1 .-<"" ' j 

_i ■ . , u , L 



,(1) 



where Oj and a e are the ion-sphere and electron-sphere 
radii, fee is the Boltzmann constant, and T is the tem- 
perature. Note, that the electron charge within an ion 
sphere exactly compensates the ion charge Zje\ Tj gives 
the ratio of characteristic electrostatic energy Z 2 e 2 /a,j to 
the thermal energy k^T. 

If Tj <C 1, the ions constitute an almost ideal Boltz- 
man gas, while for Tj > 1 they are strongly coupled by 
Coulomb forces and constitute either a Coulomb liquid 
or solid. The transformation from the gas to the liquid 
at Tj ~ 1 is smooth, without any phase transition. Ac- 
cording to highly accurate Monte Carlo calculations^ a 
classical OCP of ions solidifies at T w 175 (see, e.g. US 

We will discuss fusion reactions 
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FIG. 1: (Color online) Temperature-density diagram for 
12 C+ 16 mixture with equal number densities of 12 C and 
16 O. The dash-dot lines correspond (from top to bottom) to 
Too, Too, Fee = 1- For temperatures below these lines, the 
effects of plasma screening on corresponding reaction rates 
become strong. The dotted line is £y = 1. This line is the 
same for all i and j because of the same A/Z = 2 ratio. The 
reacting nuclei become bound along the dashed straight lines 
(the Gamow-peak energy passes through zero, see Sec. W\ . 
For lower temperatures (gray region) the mean field model 
(used in this paper) becomes inadequate. The shaded regions 
are bound by solid and long-dashed lines of constant burning 
times (1 yr and 10 6 yr). The thick dash line is the carbon 
ignition curve (Sec. IVl|l . 
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where A^° mp 
a compound nucleus 
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The reaction rate is determined 



Zi + Z-j refer to Sec. |V] and Fig. Sij is Kronecker delta, that ex- 



by an astrophysical S- factor. For a non- resonant reac- 
tion, it is a slowly varying function of energy (see, e.g., 
Ref. It is determined by the short-range nuclear 

interaction of fusing nuclei and by the Coulomb barrier 
penetration problem. The latter task can be reduced to 
the calculation of the contact probability <?r/(0), which 
is the value of the quantum-mechanical pair correlation 
function of reacting nuclei i and j at small interionic dis- 
tances r — > 0. Finally, the reaction rate can be written 
as (see Sec. [V] and (2fl|) 
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n i n S lj 5*7(0), (3) 



where 



Sij(E) = a tj (E)Eexp(2inn3) 



(4) 



is the astrophysical factor. It should be taken at an 
appropriate energy E = E^ , which is the "Gamow 
peak energy" , corrected for the plasma screening (see 



eludes double counting of the same collisions in reac- 
tions of identical nuclei (i — j). For these reactions, 
afj = h 2 / (2p,ij Zi Zj e 2 ) reduces to the ion Bohr radius. 
Furthermore, pij = m u AiAj / A c ° mv is the reduced mass 
of the nuclei and cry (E) is the fusion cross section. Fi- 
nally, rjij = ZiZje 2 yj pij/(2Eh 2 ) determines the pene- 
trability of the Coulomb barrier. 

Figure [T] is the temperature-density diagram for a 
12 C+ 16 mixture with equal number densities of both 
species. Also, we present the lines of constant dimension- 
less parameters (see Sec. Ill Bp Tij = 1 (dash-dot lines) 
and Qj = 1 (the dotted line). The dash-dot lines (from 
top to bottom) correspond to Too — lj Too = 1, and 
Tqq — !■ The Qj — 1 line is the same for all combina- 
tions of ions because of equal A/Z = 2 rations for 12 C 
and 16 O ions. The shaded regions are bounded by the 
solid and long-dashed lines of constant burning times of 
1 and 10 6 yr, respectively. The region for the O+O re- 
action is highest, for C+O intermediate, and for C+C 
lowest (as regulated by the height of the Coulomb bar- 
rier). The lines of constant burning time are similar. At 
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P IS 10 9 gem -3 they mainly depend on temperature be- 
cause of exponentially strong temperature dependence of 
the reaction rates. These parts of the lines correspond 
to thermonuclear burning (with weak or strong plasma 
screening). With growing density, the curves bend and 
become nearly vertical (that corresponds to the pyenonu- 
clear burning due to zero-point vibrations of atomic nu- 
clei; see [14j, [l(| |2lj, 122J ) . The pycnonuclear reaction rate 
is a rapidly varying function of density, which is either 
slightly dependent or almost independent of temperature. 

The carbon ignition curve (discussed in Sec. IVip is 
shown by the thick dashed line. 



B. Dimensionless parameters 

Let us consider a binary mixture of ions with mass 
numbers A% and A 2 , charge numbers Z\ and Z 2 , and 
fractional numbers X\ = n\jri\ and x 2 = 1 — x\. Follow- 
ing Itoh ct al. [23], in addition to the ion radii = an 
we introduce an effective radius 



v l/3 . ^1/3 

di + a 2 Z\ + ZJ 
ai2 = = a„ 



(5) 



2 2 
It is also convenient to define the parameter 

T c = e 2 /(a c k B T), (6) 

The coupling parameters for ions can be written as 



t a r c z^ 



5/3 



(7) 



In addition, we introduce an effective Coulomb coupling 
parameter 



Ti2 = ZiZ 2 e 2 /{a 12 k B T), 



(8) 



to be used later for studying a reaction of the type 
(A 1 ,Z 1 ) + (A 2 ,Z 2 ). 

Furthermore, we introduce the corresponding ion ra- 
dius and the Coulomb coupling parameter for a com- 
pound nucleus (A ( + Aj, Z t + Z-) = (A™ mp , Z™ mp ), 



COinp / r7COinp\ 1/3 r COmp ti / ry 

% = °e {z t3 v ) , iv = r e (z, 
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The Coulomb barrier penetrability in the classical ther- 
monuclear regime (neglecting plasma screening) is char- 
acterized by the parameter 
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The Gamow peak energy in this regime is 



(10) 



Ef = k B T nj /3. (11) 

The effects of plasma screening on the Gamow peak en- 
ergy and tunneling range are discussed in Sec. fVl 



Also, we introduce a dimensionless parameter 

Cy = 3 /rij . (12) 

In an OCP, reduces to the parameter £ defined in 
Ref. [1 71 ] and has a simple meaning - it is the ratio of the 
tunneling length of nuclei with the Gamow peak energy 
(calculated neglecting of the plasma screening effects) to 
the ion sphere radius a^-. 

III. MEAN FIELD POTENTIAL 

Direct calculation of the contact probability gij(0) is 
very complicated even in OCP. In principle, it can be 
done by Path Integral Monte Carlo (PIMC) simulations, 
but it requires very large computer resources. PIMC cal- 
culations performed so far are limited to a small num- 
ber of plasma ions (current OCP results are described in 
(2^, [25, 26]). In BIMs one needs a larger number of ions, 
Ni and N 2 , to get good accuracy. 

We will use a simple mean-field model, summarized in 
this section. It can be treated as a first approximation 
to the PIMC approach [27| . Previously we showed [TtJ it 
is in good agreement with PIMC calculations of Militzer 
and Pollock g| for OCP. 

A classical pair correlation function, calculated by the 
classical Monte Carlo (MC) method, can be presented as 



exp 



(13) 



where ify (r) is the mean field potential to be determined. 

The Taylor expansion of H%j (r) in terms of r for a 
strongly coupled ion system should contain only even 
powers [Hj]. The first two terms can be derived ana- 
lytically in the same technique as proposed by Jancovici 
[291 ] for OCP. The mean field potential at low r <C ay 
takes form 130] 



Hij(r) 



n-B^ ltr ij ^ comp 



where 



fo(Ti) + fo(Tj 



f /-pComp\ 



(14) 



(15) 



and /o(r) is the Coulomb free energy per ion in OCP. 
Here, the linear mixing rule is assumed. In principle, 
hjj can be extrapolated from MC Hij(r) data (e.g., Ref. 
[31j]). but the problem is delicate [32|] and we expect that 
using linear mixing rule is preferable. The linear mix- 
ing rule is confirmed with high accuracy [3^, [13] but its 
uncertainties were a subject of debates; see Ref. [35| for 
recent results. 

To determine the mean field potential in BIMs we per- 
formed extensive MC calculations of Hij(r). We did 129 
MC runs. In each run we calculated pair correlation func- 
tions <7n, <7i2 and g 22 for a given set of the parameters 
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FIG. 2: (Color online) Every BIM MC run is shown on both 
panels. Up: X2 versus IV Circles, pluses, crosses and dots 
refer to runs with different ratios of Z2/Z1: 1 < Z2/Z1 < 2, 
2 < Z2/Z1 < 3, 3 < Z2/Z1 < 5, 5 < Z2/Z1 < 8, respectively. 
Down: Z2/Z1 versus Vi. Circles, pluses, crosses and dots refer 
to runs with different fractions X2- X2 < 0.05, 0.05 < X2 < 0.2, 
0.2 <X2 < 0.4, 0.4 < Z2/Z1 < 1. 



Zi, Z2, N\, N2, and IV In Fig. [2] each set is shown 
by a symbol. The top panel demonstrates a fraction of 
particles 2, X2, versus IV Circles, pluses, crosses and 
dots refer to different ratios of Z2/Z1: 1 < Z2/Z1 < 2, 
2 < Z 2 /Z x < 3, 3 < Z 2 /Z x < 5, 5 < Z 2 /Z 1 < 8, re- 
spectively. In the bottom panel we show the Z2/Z1-T&H0 
versus IV Circles, pluses, crosses and dots correspond 
to different values of X2- X2 < 0.05, 0.05 < X2 < 0.2, 
0.2 < X2 < 0.4, 0.4 < X2 < 1. For some parameter sets 
(e.g., for Ti = 180, x 2 ~ 0.05, Z2/Z1 = 4/3) we use a 
number of MC runs with different initial configurations 
(random or regular lattice) and/or different simulation 
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FIG. 3: (Color online) Mean field potential uy versus r/a c 
for four MC runs. Solid lines are extracted from MC data, 
dashed lines are given by the fit (|16[) . On each plot the lines, 
which are upper at r > 2a E , correspond to U22(r); lower lines 
are 1112 and lowest lines give tin. Thin vertical dotted lines 
connected by thick solid lines show the ranges of r, which are 
used in the fitting. 
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times. We did not find any significant dependence of the 
mean field potential (for r < 2<2jj ) on the initial parame- 
ters. From all our simulations we extract -ffn(r), Hi2(r) 
and H22 (r) and remove some low-r points with large MC 
noise. We then fit these functions at r < 2a.ij by a simple 
universal analytical expression 



H l3 (r) = k B TT°f 



(16) 



Ul + U 1 x i + U 2 x s 



10 ' 



where 



y 1 + U 3 x + U 4 x 2 + U 5 x 4 + U 6 x 6 + U 2 x 
= r/aff, U Q = h%/Tf, U x = 0.0319, 



U 2 = 0.0024, 

cy— 1/3 / r-fP 1 comD\2 



u. 



-0.7 



= 0.22 (r|f) Ui 

(a?f/a-° mp ) 2 /U , U 5 = 0.111, and U 6 = 0.0305. 
Note, that Uq and U4 are fixed by the low-r asymptote 
of Hij(r) at strong Coulomb coupling, 3> 1, Eq. (fT4|) . 
For a better determination of H\2{r), we introduce the 
effective ion sphere radius af? and the Coulomb coupling 
parameter T°f: 



,cff 



ncfF 



Z l Z ] e 2 



(17) 



Here, p is an additional fit parameter which is found to 
be p = 1.6. Note, that a° = an — at for reactions with 
identical nuclei (i = j). If, however, i ^ j the parameters 
a?j and differ by a few percent and this difference is 
very important. Using to define x in Eq. (I16p and 
using Hij(r) to reproduce the pair correlation function 
gij(r) through Eq. |[T3)). we find the first peak of to 
be unrealistically deformed for large Z<zjZ\ > 5. 

Let us mention that in the OCP case our fit expression 
(fT!))) is not the same as we used previously [13] • Equa- 
tion (fT6"|) seems to be better because it better describes 
the first peak of g(r) = gu(r), especially at very large 
r > 200. For an OCP at T < 200, the old and new 
fit expressions give approximately the same accuracy of 
H(r)=H u (r). 

Four examples of MC runs are shown in Fig. [3J For 
each run (on each plot) we show dimensionless mean field 
potential Uy = Hij(r)/ (ksTTij) versus r/a e . Solid lines 
are MC data, dashed lines are derived from the fit ex- 
pression (fTB"|) . One can see numerical noise of MC data 
at low r; it is especially strong in the lowest plot for 
r o = 5. During fitting procedure we removed some low-r 
points with high MC noise. Our fitting ranges are shown 
by pairs of thin dotted lines connected by a thick solid 
line. The low-r bound of the fitting region is different 
for -Hn, H12, and H22 and for each MC run. The high-r 
bound was taken to be 2ay for all data. Thus, all our 
fits include the vicinity of the first peak of a pair cor- 
relation function at r ~ 1.8ajj. Although we have no 
MC data at low r, we expect that our approximation 
is well established because it satisfies the accurate low- 
r asymptote (fT4|). Also, it satisfies the large-r asymp- 
tote Uij(r —* 1) = aij/r, which corresponds to the fully 



screened Coulomb potential. Both asymptotes for U22 
are shown on the upper plot by dash-dot lines. One can 
see, that they nicely constrain the mean field potential. 

A poor MC statistics at low r can be, in principle, in- 
creased using specific MC schemes [36[ but this is beyond 
the scope of the present paper. 

Previously the mean field approximations was used by 
Itoh et al. [23j and Ogata et al. [3jJ. While we use the 
correct form of the linear mixing rule (fT5|) to calculate 
Mi 2 (0), the cited authors restricted themselves by OCP 
results, which depend only on r 12 (but not on Z-xjZ-y). 
The results of this simplification are discussed in Sec. lIVl 
they are visible in Fig.|4j In addition, our approximation 
is based on a more representative set of MC runs (see Fig. 
[3]) , and each run was done with better accuracy. For ex- 
ample, where is no visible noise gij (r)-noise in the vicinity 
of the first peak in out data, whereas this noise is obvi- 
ous in Fig. 1 of Ref. [37| ■ We have also checked that our 
approximation (fTo) better reproduces all our MC data 
than approximations of both groups [H, [13] • As noted 
by Ogata et al. [13], the height of the first (?y(r)-peak 
depends on the fraction of ions with larger charge (the 
peaks become higher with increasing this fraction). Our 
data confirm this statement for a not very large Coulomb 
coupling parameter, Ti < 10. Ogata et al. [37| added a 
correction for this feature in their approximation of u\\. 
It helps to describe the first gn(r)-peak for some of our 
runs (e.g., for Ti = 10, Z<xjZ\ = 5 and x 2 = 0.2). How- 
ever, we think that we still do not have enough data to 
quantitatively describe this correction to the linear mix- 
ing with good accuracy, and we use the mean field po- 
tential (|16p based on the linear mixing rule throughout 
the paper. In the Appendix we describe the corrections 
to the linear mixing on a phenomenological level. 



IV. ENHANCEMENT OF NUCLEAR 
REACTION RATES 



Let us introduce the enhancement factor F?f r of a nu- 
clear reaction rate under the effect of plasma screen- 



ing, 



ij ij ij ' 



Fl? = exp(^). 



(18) 



Here, iiy is the reaction rate in the absence of the screen- 
ing. In can be calculated using the classical theory of 
thermonuclear burning, 



nth 4 ^"j 



sum m 



3/ij 



■ exp (-nj) . 



(19) 



The barrier penetrability parameter and the Gamow 
peak energy £"? k , which enter this equation, are defined, 
respectively, by Eqs. ([TO]) and (fTT|) . 

Using Hij (r) we can calculate the reaction rate. This 
calculation takes into account the plasma screening en- 
hancement in the mean-field approximation (as was done 
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FIG. 4: (Color online) Top: Enhancement factors hij/T\2 for several Fi2 (=1, 10, 50, 200; each shown by a line on every panel) 
and Z%jZ\ (=1, 2, 5, 10; each on its own panel). The long dash-dot line presents the fit expression by Itoh et al. [S^l (IKM90); 
the short dash-dot lines are the results by Ogata et al. [33] (OIV93), calculated for Fi2 = 50, as an example (the latter lines 
are shown only for £12 < 2 and ZijZ\ < 5 which are the fit bounds in OIV93). Bottom: Ratios of R12, given by Eq. (I23|l . to 
the mean field reaction rate i?i2 F , given by l|20[l. versus £12. See text for details. 



for OCP in our previous paper 



R 



MF 




1 + 5ij 

X) 

dE exp 



knT 



(20) 



Here, E is the center-of-mass energy of the colliding nu- 
clei (with a minimum value -E m i n at the bottom of the po- 
tential well), exp(— E/k-eT) comes from the Maxwellian 
energy distribution of the nuclei, and 



Pij(E) = 



2 v /27^ 



ZiZje 2 



H i3 {r)-E (21) 



is the Coulomb barrier penetrability (r t being a classical 
turning point). 

In the WKB mean-field approximation the enhance- 
ment factor of the nuclear reaction rate is 



- R^ F {H tj } I R^ F '{Hij - 0}. 



(22) 



In the cases of physical interest, R^ {0} can be inte- 
grated by the saddle-point method. However, we always 
integrate in Eq. (f2"0"|) numerically. 



The calculated enhancement factor Fff r can be fitted 



as 



logF, 



scr 

ij 



+ /o 



i-icomp 



(23) 



where /o(r) is a free energy per ion in an OCP. We use the 
most accurate available analytic fit for /o(r) suggested 
by Potekhin and Chabrier [381 ] : 



/ (r) = a 1 Vr (A 2 + t) 



A 2 In 

2A 3 

Bi 
B 3 



— arctan Vr 



In 1 



-02 In ( 1 

11 

B 4 



r 

~~Eh, 



(24) 



= -0.907, A 2 = 0.62954, A 3 = -a/3/2 - 
i 0.2771, B x = 0.00456, B 2 = 211.6, B 3 = 
— 10~ 4 , B 4 = 0.00462. We also introduce a number of 
additional parameters: 



Here, Ai 



4 Zi Z* 



('2 



{Zi + z 

0.406 y. 



ci = 0.013 y. 



3) 

0.14 
ij -i 



C3 



0.062 1/9: 19 + i.8/r 



tij = (1 • ' • <:4) 



1/3 



(25) 



The root-mean-squared relative error of the fit of F scr 
is 5%. The maximum error ~ 14% occurs at Qj w 4.8, 
T 12 w 25, and Z 2 /Zi sa 10. The fit was constructed for 
1 < Z 2 jZ x < 10, 1 < r 12 < 200 (we also use r 12 = 400 
and 600 for testing) and £12 < 8 (larger values of C12 were 
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not included into fitting but used for tests). As for OCP 
in our previous paper [17j . at larger T the relative errors 
increase. 

In the limit of Cu -> we have log_K?f r = h%. This is 
the well known relation (e.g., [3JJ |40J ) ; its formal proof 
is given in [41| . 

Our fit formula (|23[) has the same form as the thermo- 
dynamic relation (|15[) . but contains Coulomb coupling 
parameters divided by Ty . It generalizes our OCP 
result (l7j . 

Figure 3] shows the normalized enhancement factor 
h\2/T\2 (bottom panels) as a function of £12 for several 
values of Ti2 (= 1, 10, 50, 200; separate lines on each 
plot) and Z2/Z1 (= 1, 2, 5 and 10, separate plots). The 
first panel (for Z2 = Z\) corresponds not only to a re- 
action in OCP, but to reactions of identical nuclei in a 
mixture; the mean field potential and the plasma screen- 
ing enhancement do not depend on an admixture of other 
elements [see Eqs. (TIBl) and (|2"5j) ]. as long as the Coulomb 
coupling parameter is not too small, > 1, and the lin- 
ear mixing is valid. The larger Ty , the weaker is the de- 
pendence of hyxfVyx on Tij. Long-dashed lines, plotted 
for Ti2 = 50, significantly differ from dashed lines which 
are for Ti2 = 10, but can be hardly distinguished from 
solid lines, plotted for Ti2 = 200. Therefore, the total en- 
hancement factor exp (/112) depends exponentially on Ti2 
in the first approximation. The normalized enhancement 
factor hi 2 /Ti2 suggested by Itoh et al. [H| (IKM90) is 
shown by long dash-dot lines. It is independent of T12, 
and we show one line in each plot. The short dash-dot 
lines demonstrate the fit expression of Ogata et al. [37| 
(OIV93) (for Ti2 = 50 as an example). The lines are cut 
at C12 = 2 (that bounds the fit validity). We show no 
OIV93 line on the panel for Z2/Z1 = 10 because there 
are no OIV93 simulations for such large Z2/Z1 ratios. 
The two lines, IKM90 and OIV93, demonstrate higher 
enhancement, especially for (12 < 1 and large Z2/Z1. 
First, IKM90 and OIV93 used less accurate thermody- 
namic approximations for calculating Uij(0). For large 
Z2/Z1 -ratios, their fit error increases because of using 
the OCP expression (dependent only of Fi 2 , but not on 
Z>2 / Z\ ) for determining Uij (0) . This increases a difference 
between our and their approximations at large Z2 jZ\ (see 
also a discussion in Sec. IIIC of Ref. [3). IKM90 and 
our (^-dependence of the enhancement factors is quali- 
tatively the same, but the results of OIV93 (and of the 
preceding paper [3fJ) are qualitatively different. As we 
show previously |l7| . such results contradict recent PIMC 
results by Militzer and Pollock [2(| for OCP. The devia- 
tion of the IKM90 enhancement factors at low (y values 
comes possibly from using an oversimplified mean field 
potential. 

The indicated differences between IKM90, OIV93 and 
our results can lead to very large differences of the total 
enhancement factor exp(/ii2). For example, at ^2/^1 = 
5 and Ti2 = 200 the difference can reach five orders of 
magnitude. Accordingly, we do not present the IKM90 
and OIV93 results in lower panels, which give the ratio of 



the reaction rates Rf^, given by the fit expression (|23[) . 
to the reaction rate i?^ F , calculated in the mean field 
approximation using Eq. (f2"0"|) . A very small difference 
(< 10%), shown on these plots, is true only by adopting 
the mean field potential (fT6|) . The uncertainties of the 
reaction rates, which come from the uncertainties of the 
mean field potential and inaccuracy of the mean field 
approximation, can be larger. 

Let us stress that Eq. (|2"3"|) is derived for the case of 
strong screening; it is invalid at Ti2 < 1, where the 
adopted linear mixing rule fails. A simple correction of 
the enhancement factor for weaker screening is suggested 
in the Appendix. It allows us to extend the results for 
the case of weaker Coulomb coupling. 



V. ASTROPHYSICAL S-FACTORS: GAMOW 
PEAK ENERGIES AT STRONG SCREENING 

In the previous section we focused on the plasma 
screening of the Coulomb barrier penetration in dense 
plasma environment. Now we discuss the effects of 
plasma screening on the astrophysical S'-factor, that de- 
scribes nuclear interaction of the reacting nuclei after the 
Coulomb barrier penetration. The latter effects are not 
expected to be very strong because of not too strong en- 
ergy dependence of Sij(E) (e.g., Ref. [13 )> but we con- 
sider them for completeness. 

The astrophysical S'-factor describes the effects of 
short-range nuclear forces and it should depend on the 
parameters of the reacting nuclei just before a reaction 

event. The energy Ef^' substituted into S(E) in Eq. pO]) 
should be corrected for the mean field potential created 
by other ions, £?■ = E + Hij (0). This correction is obvi- 
ous and has been used in calculations (e.g., [HI]), but its 
formal proof has not been published, to the best of our 
knowledge. The reaction rate in the presence of plasma 
screening can be calculated as 

mnj ( 8 \ 1/2 
" ~ 1 + % WijT^J 

x Ja ij (E)Eexp(-^jdE, (26) 

where <7y (E) is the reaction cross section including the 
screening effects. To calculate Rij, let us remind the bar- 
rier penetration model a with parameter-free model of 
nuclear interaction gives a good description of reaction 
cross sections [l^. In the absence of plasma screening 
at not too high energies E (where only the s-wave chan- 
nel is important), this model reduces to the WKB cal- 
culation of the penetrability through a potential which 
is the sum of the Coulomb potential and the potential 
Vy(r,v?j), that describes the short-range nuclear inter- 
action. The latter potential depends on the local relative 
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velocity v(r.E) of the reactants given by 
o 2 



My 



E + H l3 {r)-^^-Vmry) 



(27) 



Let us generalize this model by adding the screening po- 
tential Hij (r) and write the cross section as 



7T 



2^/2/iy 



dr 



(28) 



Here, r tn and r t are classical turning points and k 2 j = 
2/dijE/H 2 . Substituting this equation into Eq. (|26|) we 
have 



1? 



1 + 5,j \TTfJ, i3 T 3 



1/2 



x / Sij(E)exp 



E 
T 



p ij( E ) 



dE, (29) 



where Sij(E) is the astrophysical factor, calculated in the 
presence of the plasma screening, 



S ij (E)=a ij Eexp[P ij {E)}; 



(30) 



Sij(E) — Sij(E + Hij(Q)) and the reaction rate is 



R 



1/2 



1 + Sij \nfi i: jk^T 3 
J S l3 (E + H l3 (0))e X p 



E 



d^32) 



If this simple barrier penetration theory were invalid 
and the hypothesis of low-energy hindrance of nuclear 
reactions [44| were correct, such a simple correction of 
astrophysical S-factors for the plasma screening effects 
could be insufficient, and the calculation of the reaction 
rate would be much more complicated. We will not con- 
sider this possibility. 

As in of the absence of the plasma screening, the main 
contribution to the integral (|32[) comes from a narrow 
energy range. We neglect the energy dependence of 
Sij(E + Hij(0)) in this range and take the S'-factor out 
of the integral. The contact probability in the mean field 
model is given by 



2^FJI~Z 1 Z 2 e 2 

9iM= hTt^ » / " xi> 



E 
k^T 



dE, 
(33) 



so that the reaction rate reads 

+*«<»>)*<<>>• < M > 



P l3 {E) is defined by Eq. (|2T1)- The potential V$(r,v) is 
non-zero only for very small r < rtj. Let us substitute 
Eqs. dHJ) and fl55j! into flH|) and divide all integrals into 
two parts, at r < and r > r^. The integrals over r > 
tn which come from Pij(E) and aij{E) will be exactly 
the same (V^(r, v) = for r > r^) and cancel each other. 
As a result, 



Sij (E) = 7r — cxp t 



E 



We have calculated the modified Gamow-peak energy 
? k = Ef^ + Hij(0), which should be used as an argu- 



ment of S'-factors, in a wide range of plasma parameters. 
It can be fitted as 



E 



P k' 



Hi 



(0) 



k B T^ 



1/3 



(35) 



2/ij 



ZiZje 2 



Hij{r)-Edr 



(3i; 



The maximum relative fit error for 1 < L.y < 200, Qj < 8 
and 0.1 < Zi/Zj < 10 is 4%. It takes place at T l3 = 1, 
Qij = 1 and Zi = Zj. Because of large nuclear physics 
uncertainties in our knowledge of S(E) at low energies of 
astrophysical interest, this accuracy is more than suffi- 
cient. An example of the dependence of E± 2 /^12 on C12 
is shown in Fig.[Ufor Zq/Zj. = 2 and T 12 = 100. The solid 
line is the fit (|35|) . the dashed line is a result of the ex- 

act £fy calculation in the mean field model. The peak 

energy E^- has two asymptotes shown by thin dotted 
V£j (r, v 2 ) in the low-r region r < r N , where the nuclear lines. At low C < 1, the reaction occurs in the thermonu- 



fl£jf_ _ Hi j (r) + V$(r, v 2 )-Edr 



As expected, the astrophysical S'-factor is determined by 
the behavior of the total potential ZiZjt 2 jr — Hj(r) + 



forces are important. Because tn is much smaller than 
the ion sphere radius ay , which is a typical length scale of 
Hij(r) (e.g., Fig. [3]), we can neglect variations of Hj{r) 
for r < tn and replace i/y(r) under the integrals in 
Eq. by .Hy(O). The Sij (J5)-factor in the absence 
of plasma screening is given by the same Eq. (|31[) but at 
Hij(r) = 0. Thus, Sij(E) is defined by the same equation 
as Sij(E'), provided E' = E + ffy(0). In other words, 



clear regime, Ef^ 

sical result. For large Cij 3> 1, we have E, 



^ij/Cij: which is the standard clas- 



Hij(0). 



The classical asymptote Ef^ = ry/£y is applicable at 
dj 0-5, where the plasma screening enhancement can 
be very strong (tenths orders of magnitude). This is be- 
cause in the thermonuclear regime the tunneling length is 
not very large and Hij (r) does not significantly change it 
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FIG. 5: (Color online) The Gamow peak energy Ef^ normal- 
ized for Ti2 as function of £12 at Z2/Z1 = 2 and Fi2 = 100. 
Solid line is the fit expression (|35[l . dashed line shows the 
mean field result. Thin dotted lines correspond to the low-^j 

thermonuclear asymptote Ef^ = Ti2/Ci2 and high-^j pyc- 

nonuclear asymptote Ef^ = Hij(0). See text for details. 



(see Fig. [6J; ions tunnel through the Coulomb potential, 
shifted by i?y(0). The shift increases the probability of 
close ions collisions by a factor of exp (if^ (0)/T), but 
does not change the Maxwellian energy distribution and 
the dependence of the tunneling probability on tunneling 
length. As a result, the main contribution to the reaction 
rate comes from ions with the same tunneling length (and 

the same Ef^ ), as in the absence of plasma screening. 

The large-(y asymptote is simple. For such parame- 
ters, thermal effects are small and the ions, which mainly 
contribute to the reaction rate, correspond to the min- 
imum energy of the total potential Z\Zijr + Hij(r). 
This energy is small compared with ffy(0), so that 

Ef^ w Hij (0) . One can see a small difference of Hij (0) , 
shown by the thin dotted line in Fig. O and the asymp- 
tote of Hij(Q) at large This difference does not in- 
troduce significant uncertainties to the reaction rate be- 
cause of much larger nuclear-physics uncertainties caused 
by our poor knowledge of S(E) at low energies (e.g., Ref. 

(TBI). We assume that the Ef^ approximation ([55]) is 
valid (at least qualitatively) not only for thermonuclear 
burning with strong screening, but also for pycnonuclear 
burning. In the pycnonyclear regime, the ions occupy 
their ground states and oscillate near their lattice sites. 
The reaction rates are determined by zero-point vibra- 
tions of the ions. The energy of zero-point vibrations 
is typically small compared to the minimum energy of 



the potential ZiZj/r + Hij(r) (the latter is much lower 
than Hij(0), where Hij{r) is an anisotropic effective po- 
tential created by neighboring ions). Therefore, the ions 
start with a small kinetic energy (the minimum of the 
ZiZj/r + Hij(r)) and tunnel to r — + through the 
Coulomb potential plus the Hij(r). During tunneling, 
they fall into the potential well Hij(r), and their energy 
increases by Hij(Q). Therefore, Ef^ should be equal to 
Hij(0), but not to energy of zero-point vibrations, as as- 



Note, that in the relaxed-lattice 
= fly(0) 



sumed in 01, [13 
approximation ily(O) 

For example, let us consider Ef^ and the characteris- 
tic (half-maximum) energy widths of the Gamow peak 
for a 12 C and 16 O mixture (with equal number densities 
of C and O ions) at p — 5 x 10 9 g cm~ 3 . In Fig. [S] we 
show the effective total radial mean-field Coulomb po- 



tentials Uijir) 



.12 



C or 



16 0) 

for five temperatures, T = 10 9 , 10 85 , 10 8 , 10 75 , and 
10 7 K. For this mixture, the ion-sphere radius of 12 C 
ions is ctc=98 fm, and it is ao=108 fm for 16 ions. The 
solid, dotted and long-dash lines show the CC, CO and 
OO potentials, respectively. Each potential Uij{r) has 
a minimum at 2atj due to the Coulomb coupling. 

The thin horizontal lines connected by double-arrow lines 
in Fig.[T5]show the Gamow-peak energy ranges. The types 
of thin lines are the same as for [Zy(r). Thick sections of 
short horizontal lines, which intersect double-arrow lines, 
demonstrate the Gamow-peak energies Ef*. Thick parts 
of short-dashed lines at small r < 100 fm indicate the 
thermal energy level k^T measured from the bottom of 

The top right panel marked "without screening" shows 
the same lines (as other panels) for T = 10 s K neglect- 
ing plasma screening. Dramatic difference of unscreened 
and screened potentials and corresponding Gamow-peak 
regions is obvious from comparison with bottom panels. 
The plasma screening reduces the Gamow-peak energy 
Ef^ by a factor of four, but (as described above) does 
not affect significantly the Gamow-peak width. 

Let us discuss Fig. |J5] in more detail. The panels plot- 
ted for T = 10 8 - 5 and 10 9 K (top left and top center, 
respectively) refer to the thermonuclear reaction regime 
with strong plasma screening. The next two panels 
(right bottom and center bottom) are for a colder plasma 
(T = 10 8 K and T = 10 7 5 K, respectively), while the last 
(left bottom) panel is for a very cold plasma (T = 10 7 K) 
(that is certainly in the zero-temperature pycnonuclear 
regime). When the temperature decreases, the Gamow- 
peak energy range becomes thinner (note the difference 
of energy scales in different panels) and shrinks to lower 
energies. If T > 10 8 K, the Gamow peak range is still 
at E > [belonging to continuum states in a potential 
Uij (r)] , and the peak energy Ef^ (the short-dash horizon- 
tal line) is in the center of the Gamow-peak range. The 
energies within this range are much higher than k^T sup- 
porting the statement that the main contribution to re- 
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FIG. 6: (Color online) Effective mean-field Coulomb potentials E/y(r) for the reactions 12 C + 12 C (solid line), 12 C + 1G O 
(dotted line) and 16 + 16 O (long-dash line) in the mixture of 12 C and 16 with equal number densities of these nuclei at 
p = 5 x 10 9 g cm -3 and five temperatures (log 10 T[K] = 9, 8.5, 8, 7.5, and 7). 



action rates at sufficiently high T comes from suprather- 
mal ions. In these cases, the underlying mean-field WKB 
approximation is expected to be adequate. In the forth 
panel, the lowest energies of the Gamow-peak range be- 
come negative (drop to bound states) and Ef^ moves 
from the center to the lower bound of the Gamow-peak 
range (so that the Gamow peak becomes significantly 
asymmetric). The mean-field WKB approach based on 
the spherically symmetric mean field potential ffy(r) 
may be still qualitatively correct but becomes quantita- 
tively inaccurate. Note that the Uij (r)-potential becomes 
positive for r > 200 fm. This feature is determined by 
the right wing of the first <7y(r)-peak, where gij{r) be- 
comes smaller than one. This wing is unimportant in 
the mean field model; our approximation gives qualita- 
tively correct Uij(r) for r > 2a%j {gij{r) — > 1 for r — > oo). 
For the lowest temperature in Fig. [5] the Gamow-peak 
energy range fully shrinks to bound-state energies and 
the formal Gamow-peak energy becomes lower than ksT, 
nearly reaching the lower bound of the Gamow-peak re- 
gion. The mean-field WKB approximation breaks down 
at these low temperatures, and the formally calculated 
£"? k is inaccurate. Nevertheless, the energy in the argu- 
ment of the S-factor, Eff' = £? k + i?y(0) » i?y(0), is 
expected to be well defined. 



VI. RESULTS AND DISCUSSION 



Figure [7] shows the temperature dependence of the 
three (C+C, C+O, and O+O) reaction rates in a 12 C 
and 16 O mixture with equal number density of carbon 
and oxygen ions at p = 5 x 10 9 g cm~ 3 . The solid 
lines are our mean-field calculations (marked as MF), the 
dashed lines are given by the fit expression (|2"5|) (marked 
as Fit), and the dotted lines are calculated neglecting the 
plasma screening (_ffjj(r) = 0). Lines of the same type 
refer (from top to bottom) to the C+C, C+O, and O+O 
reactions, respectively. The right vertical scale gives typ- 
ical carbon burning time tqc — nc/Rcc- The shaded 
region (T < 2.5 x 10 7 K) corresponds to bound Gamow- 
peak states (see Fig. [5] and corresponding discussion in 
the text). This region is similar for all three reactions 
because of approximately the same charges of reacting 
ions. For such low temperatures the mean field model 
with isotropic potential is quantitatively invalid but qual- 
itatively correct; the reaction rates become temperature 
independent, which is the main property of pycnonuclear 
burning. The pycnonuclear burning rates are rather un- 
certain [lfij : typical rates reported in the literature are 
one order of magnitude smaller than those extracted from 
our mean field calculations at T — > 0. The difference 
may result from the fact that we use spherically sym- 
metric mean-field potential (rather than more realistic 
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FIG. 7: (Color online) Reaction rates in a 12 C+ 16 mixture 
with xc = xo — 0.5 at p — 5 x 10 9 g cm -3 . Solid lines 
are the mean-field calculations, dashed lines correspond to 
our fit expression (|23[l , and dotted lines are obtained neglect- 
ing the plasma screening. The lines of the same type refer 
(from top to bottom) to the C+C, C+O, and O+O reactions. 
The shaded region of T < 2.5 x 10 7 K corresponds to bound 
Gamow-peak states, Ef^ < 0. The right vertical scale shows 
characteristic carbon burning time rcc = nc/Rcc- See text 
for details. 



anisotropic potential). In reality, ions can be localized in 
deeper potential wells near their lattice sites. 

The fit and mean-field results, which are mostly in- 
distinguishable in the figure, can strongly differ in the 
shaded region, especially for the O+O reaction. We do 
not expect that it is a significant disadvantage of our 
fit expression (because the mean field approach becomes 
invalid at such conditions), but we would like to men- 
tion this feature. Note, that for the O+O reaction in 
Fig. [7] the temperature T — 2.5 x 10 7 K corresponds to 
Too ~ 400 and the total enhancement factor is ~ 10 110 . 

Figure [5] presents carbon ignition curves, which are 
most important for modeling nuclear explosions of mas- 
sive white dwarfs (supernova la events) and carbon explo- 
sions in accreting neutron stars (superbursts). For white 
dwarfs it is determined as a line in the T — p plane, where 
the nuclear energy generation rate equals the local neu- 
trino energy losses (which cool the matter). For higher 
T and p (above the curve), the nuclear energy genera- 
tion exceeds the neutrino losses and carbon ignites. The 
curves are plotted for 12 C+ 16 mixtures. The main en- 
ergy is generated in the C+C reaction even for xc — 0.01 
because the C+O and O+O reactions are stronger sup- 
pressed by the Coulomb barrier (see Fig.[7]to compare the 
reaction rates). In our mean-field model, which employs 
the linear mixing, the admixture of oxygen affects the 
C+C burning only by reducing the number density of car- 
bon nuclei [not through the contact probability gcc(0); 
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FIG. 8: (Color online) Carbon ignition curve in 12 C+ ie O 
mixtures. Dash-dotted, solid and dashed lines show the mean- 
field calculations for xc = 1, 0.5, and 0.01, respectively. We 
also plot two thin dotted lines of constant rcc = 1 and £cc = 
1 to demonstrate typical values of dimensionless parameters. 
See text for details. 

see Eq. (|34p]. so that the C+C reaction rate is oc x c . The 
neutrino energy losses are mainly produced by plasmon 
decay and electron-nucleus bremsstrahlung processes. 
The neutrino emissivity owing to plasmon decay is cal- 
culated using the results of Ref. [45| (with the online ta- 
ble http://www.ioffe.ru/astro/NSG/plasmon/table.dat). 
The neutrino bremsstrahlung emissivity is calculated us- 
ing the formalism of Kaminker et al. (4o| . which takes 
into account electron band structure effects in crystalline 
matter. For an CO mixture, this neutrino emissivity is 
determined using the linear mixture rule. 

Two thin dotted lines in Fig. [5] correspond to constant 
Tec = 1 an d Ccc = 1- Dash-dotted, solid and dashed 
lines are the mean-field calculations for xq = 1 (pure car- 
bon matter), 0.5, and 0.01, respectively. For a fixed p, the 
higher xc the higher the number density of carbon ions 
and the higher the reaction rate. This intensifies carbon 
burning, and the carbon ignites at lower temperature. 
For high densities, the carbon ignition curve bends and 
the carbon can ignite at very low temperatures. This is 
caused by weakening the temperature dependence of the 
reaction rate (see Fig. [7]) and by a strong suppression of 
neutrino emission with decreasing temperature. 

VII. CONCLUSIONS 

We have analyzed the plasma screening enhancement 
of nuclear reaction rates in binary ionic mixtures. We 
have used a simple model for the enhancement factor 
based on the radial WKB tunneling of the reacting nu- 
clei in their Coulomb potential superimposed with the 
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needs to extrapolate experimental results to lowest ener- 
gies. Throughout the paper we assumed a smooth energy 
dependence of astrophysical S-factors, that is supported 
by calculations in the frame of the barrier penetration 
model (e.g., Ref. [42|). However, some models (e.g., 
predict resonances at low energy, which can significantly 
change the reaction rates and ignition curve (48| . We do 
not discuss such effects in the present paper. New exper- 
imental and theoretical studies of astrophysical S'-factors 
are needed to solve this problem. 
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static mean-field potential created by neighboring plasma 
ions. We have done accurate Monte Carlo calculations 
of the mean-field plasma potential for a two-component 
strongly coupled plasma of ions and proposed a simple 
and accurate analytic fit to the plasma potential (Sec. 
IIIip . We have calculated the plasma enhancement fac- 
tors of nuclear reaction rates in the mean-field WKB ap- 
proximation and have obtained their accurate fit (Sec. 
IIV[) . We have analyzed the effect of the plasma screen- 
ing on astrophysical S'-factors and Gamow-peak energies 
(Sec. [V}. To illustrate the results, we analyzed nuclear 
burning in 12 C+ 16 mixtures fSec. IVI[) . 

We demonstrate that the mean-field WKB method 
gives qualitatively correct (temperature independent) re- 
action rates even in the zero-temperature pycnonuclear 
burning regime. In this regime the dynamics of the re- 
acting ions is determined by zero-point vibrations; they 
fuse a long selected (anisotropic) close-approach trajec- 
tories [14]]; the mean- field radial WKB method was ini- 
tially expected to be absolutely inadequate. Let us men- 
tion in passing that the problem of pycnonuclear burn- 
ing has not been accurately solved even for OCP plasma 
(e.g., see [13), and the uncertainties of the solution in- 
crease in multicomponent mixtures (see (l6j . and refer- 
ences therein). 

Other uncertainties in our knowledge of the reaction 
rates come from nuclear physics. As a rule, the astro- 
physical S-factors cannot be experimentally measured 
for such low energies as Gamow peak energies in stel- 
lar matter. For example, the lowest experimental point 
for the 12 C+ 12 C reaction is ~ 2.1 MeV, whereas typical 



Gamow peak energies Ef, are ~ 1 MeV. Therefore, one where 



APPENDIX A: CORRECTION OF THE 
ENHANCEMENT FACTORS FOR WEAK 
NONIDEALITY 

Our fit expression (f2"3")l does not reproduce the 
well-known Debye-Hiickel enhancement factor for weak 
Coulomb coupling (Ty -C 1), 



h^ = i 1/2 Z l Z 3 (Z 2 ) 1/2 Tl' 2 /{Z) 



a/2 



(Al) 



Let us remind that is related to the reaction rates 
by Eq. (|18p . To correct our approximation in the weak- 
coupling limit we note, that a formal use of the linear 
mixing rule [employed in Eq. (|23p ] at low T^j <C 1 gives 



-1/2 



,5/2 



7 5/2 



7 5/2 



T 3 J 2 . (A2) 



Therefore, our fit (|2"3"|) gives the correct power low (oc 

To^ 2 ), but an inexact prefactor (for 0.1 < Z\jZ% < 10 
and all X\ the relative error does not exceed 40%). We 
suggest to introduce a correction factor 



a 



(z 2 ) 



1/2 



1 



(z) 1 '* (z l + z,) 5/2 ~z\ 12 -z) 12 '' 

and finally write the enhancement factor as 




produces the correct Debye-Hiickel asymptote of the en- 
hancement factor in the weak coupling limit (Tij « 1) 
and it reproduces also our results at T > 1. 

Our interpolation expression is certainly a simplifica- 
tion, but it is expected to be qualitatively correct. In 
Fig. [9] we show the enhancement factors h%j , normalized 

3/2 

with respect to r o , as a function of r c . We take a BIM 
with Z2/Z2 and X2 = 0.05, 0.5 and 0.95 as an example. 
The thin dotted lines correspond to the Debye-Hiickel 
asymptote at T <C 1. The thin solid, dash and long-dash 
lines are the linear mixing asymptotes (h{f, h 1 "^, and 
h 1 ^ , respectively), which are valid in the limit of strong 
Coulomb coupling, T ^> 1. The thick lines show the 
enhancement factors, calculated using our interpolation 
expression (IA4p . One can see that the asymptotes fix 
the enhancement factors quite well and our interpolation 
looks reasonable. 

Of course, an accurate description of the enhancement 
factors at moderate Coulomb coupling is desirable but 
the exact solution may be complicated. Fortunately, in 
all cases of physical interest the reactions in this regime 
occur at ■< 1, so that we can calculate the enhance- 
ment factors as a difference of free energies before and 
after a reaction event. Recent calculations of Potekhin 
et al. [HI of the free energy of binary and triple mixtures 
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FIG. 9: (Color online) The enhancement coefficients hij /T^ 2 (at small fy) as a function of T c for binary mixtures with 
Z2/Z1 = 2 and X2 = 0.05 (left), 0.5 (center) and 0.95 (right). Thin dotted horizontal lines correspond to Debye-Hiickel 
enhancement (|A1|) . Thin solid, dash and long-dash lines are h\f, hS, and h\™, respectively. Thick lines demonstrate the 
enhancement factors, calculated using our interpolation expression (|A4fl . 



at intermediate Coulomb coupling seem to be most im- 
portant in this respect. However, they are not very con- 
venient for practical purpose because the enhancement 
factors depend on derivatives of the free energy with re- 



spect to numbers of particles. These derivatives should 
be calculated at small concentrations of compound nuclei, 
and the accuracy of numerical differentiation deserves a 
special study. 



D. D. Clayton, Principles of Stellar Evolution and Nucle- 
osynthesis (University of Chicago Press, Chicago, 1983). 
P. Hoflich, Nucl. Phys. A 777, 579 (2006). [16 
T. Strohmayer and L. Bildsten, in Compact Stellar X- 
Ray Sources, edited by W. H. G. Lewin, M. Van der Klis [17 
(Cambridge University Press, Cambridge, 2006), p. 113. 
A. Cumming, J. Macbeth, J. J. M. in 't Zand, and D. [18 
Page, Astrophys. J. 646, 429 (2006). 
S. Gupta, E. F. Brown, H. Schatz, P. Moeller, and K.-L. [19 
Kratz, Astrophys. J. 662, 1188 (2007) 

D. G. Yakovlev, K. P. Levenfish, Contrib. Plasma Phys. [20 
43, 390 (2003). [21 
D. G. Yakovlev, K. P. Levenfish, and P. Haensel, Astron. 
Astrophys. 407, 265 (2003). [22 
D. G. Yakovlev, K. P. Levenfish, A. Y. Potekhin, O. Y. 
Gnedin, and G. Chabrier, Astron. Astrophys. 417, 169 [23 
(2004). 

D. G. Yakovlev, K. P. Levenfish, O. Y. Gnedin, Eur. [24 
Phys. J. A 25, 669-672 (2005). [25 
D. Page, U. Geppert, and F. Weber, Nucl. Phys. A777, 
497 (2006). [26 
K. P. Levenfish and P. Haensel, Astrophys. Space Sci. 
308, 457 (2007). [27 
P. S. Shternin, D. G. Yakovlev, P. Haensel, and A. Y. 
Potekhin, Mon. Not. R. Astron. Soc. 382, L43 (2007). [28 

D. G. Yakovlev, L. Casques, M. Wiescher, Mon. Not. R. [29 
Astron. Soc. 371, 1322 (2006). [30 

E. E. Salpeter and H. M. Van Horn, Astrophys. J. 155, 

183 (1969). [31] 
L. R. Gasques, A. V. Afanasjev, E. F. Aguilera, M. 



Beard, L. C. Chamon, P. Ring, M. Wiescher, and D. G. 
Yakovlev, Phys. Rev. C 72, 025806 (2005). 

D. G. Yakovlev, L. R. Gasques, A. V. Afanasjev, M. 
Beard, and M. Wiescher, Phys. Rev. C 74, 035803 (2006). 

A. I. Chugunov, H.E. DeWitt, D.G. Yakovlev, Phys. Rev. 
D 76, 025028 (2007). 

J. W. Negele, D. Vautherin, Nucl. Phys. A 207, 298 
(1973). 

H. E. DeWitt, W. Slattery, D. Baiko, andr D. Yakovlev, 

Contrib. Plasma Phys. 41, 251 (2001). 

S. Ichimaru, Rev. Mod. Phys. 65, 255 (1993). 

H. M. Van Horn and E. E. Salpeter, Phys. Rev. 157, 751 

(1967). 

S. Schramm and S. E. Koonin, Astrophys. J. 365, 296 
(1990); 377, 343(E) (1991). 

N. Itoh, F. Kuwashima, and H. Munakata, Astrophys. J. 
362, 620 (1990). 

S. Ogata, Astrophys. J. 481, 883 (1997). 

E. L. Pollock and B. Militzer, Phys. Rev. Lett. 92, 021101 

(2004) . 

B. Militzer, E. L. Pollock, Phys. Rev. B 71, 134303 

(2005) . 

A. Alastuey and B. Jancovici, Astrophys. J. 226, 1034 
(1978). 

B. Widom, J. Chem. Phys. 39, 2808 (1963). 
B. Jancovici J. Stat. Phys. 17, 357 (1977). 

S. Ogata, H. Iyetomi, and S. Ichimaru, Astrophys. J. 372, 
259 (1991). 

H. DeWitt and W. Slattery, Contrib. Plasma Phys. 39, 
97 (1999). 



14 



[32] Y. Rosenfeld, Phys. Rev. E 53, 2000 (1996). 
[33] H. DeWitt, W. Slattery and G. Chabrier, Physica B 228, 
21 (1996). 

[34] H. DeWitt and W. Slattery, Contrib. Plasma Phys. 43, 
279 (2003). 

[35] A. Y. Potekhin, G. Chabrier, F. J. Rogers, Phys. Rev. E 

79, 016411 (2009). 
[36] J.-M. Caillol and D. Gilles, J. Phys. A 36, 6243 (2003). 
[37] S. Ogata, S. Ichimaru, and H. M. Van Horn, Astrophys. 

J. 417, 265 (1993). 
[38] A. Y. Potekhin and G. Chabrier, Phys. Rev. E 62, 8554 

(2000). 

[39] E. E. Salpeter, Aust . J. Phys. 7, 373 (1954). 

[40] H. E. DeWitt, H. C. Graboske, and M. S. Cooper, As- 
trophys. J. 181, 439 (1973). 

[41] D. G. Yakovlev and D. A. Shalybkov, Soviet Sci. Rev. 
Sec. E 7, 313 (1989). 



[42] L. R. Gasques, A. V. Afanasjev, M. Beard, J. Lubian, 
T. Neff, M. Wiescher, and D. G. Yakovlev, Phys. Rev. C 
76, 045802 (2007). 

[43] N. Itoh, N. Tomizawa, S. Wanajo, and S. Nozawa, As- 
trophys. J. 586, 1436 (2003). 

[44] C. L. Jiang, K. E. Rehm, B. B. Back, and R. V. F. 
Janssens, Phys. Rev. C 75 015803 (2007). 

[45] E. M. Kantor, M. E. Gusakov, MNRAS 381, 1702-1710 
(2007). 

[46] A. D. Kaminker, C. J. Pethick, A. Y. Potekhin, V. Thors- 
son, and D. G. Yakovlev, Astron. Astrophys. 343, 1009 
(1999). 

[47] R. Perez-Torres, T.L. Belyaeva, E. F. Aguilera, Physics 

of Atomic Nuclei 69, 1372 (2006). 
[48] R. L. Cooper, A. W. Steiner, and E. F. Brown, 

larXiv:0903.3994l 



